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Abstract 

Modifications on the predictions about the matter power spectrum based on the hypothesis of 
a tiny contribution from a degenerate Fermi gas (DFG) test-fiuid to some dominant cosmological 
scenario are investigated. Reporting about the systematic way of accounting for all the cosmolog- 
ical perturbations, through the Boltzmann equation we obtain the analytical results for density 
fluctuation, 6, and fluid velocity divergence, 9, of the DFG. Small contributions to the matter 
power spectrum are analytically obtained for the radiation-dominated background, through an 
ultra-relativistic approximation, and for the matter-dominated and A-dominated eras, through a 
non-relativistic approximation. The results can be numerically reproduced and compared with 
those of considering non-relativistic and ultra-relativistic neutrinos into the computation of the 
matter power spectrum. Lessons concerning the formation of large scale structures of a DFG 
are depicted, and consequent deviations from standard ACDM predictions for the matter power 
spectrum (with and without neutrinos) are quantified. 
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I. INTRODUCTION 



The cosmic accelerated expansion related to the conception of dark energy, the nature 
of dark matter components, the role of cosmological background neutrinos, and finally, the 
comprehension of all their intrinsic relations belong to one of the most challenging current 
problems in theoretical physics. Such an enlarged overview of the dark sector does necessarily 
involve a fruitful interplay between general relativity, astrophysics and particle physics [1-4] , 
which is fundamental to the linear perturbation theory that describes the cosmic inventory. 

Besides presenting a consistent theoretical foundation, the theory for cosmological pertur- 
bations has been highly predictive. In particular, it has explained the precise measurements 
of temperature and polarization anisotropics of the cosmic microwave background (CMB), 
most notably those from the Wilkinson Microwave Anisotropy Probe [5]. Even being a 
matter of the same scope, the fine-tuning between theoretical predictions and observable 
data for large scale structures deserves, however, a sharper analysis. In fact, it falls on more 
complex procedures for which even the nonlinearity effects are sometimes relevant [6] . The 
matter power spectrum describes the density contrast of the Universe, which parameterizes 
the inhomogeneities in the matter distribution of the Universe. The density contrast corre- 
sponds to the difference between the local energy density and the averaged energy density 
of some component of the cosmic inventory. On large scales, the gravity competes with the 
cosmic expansion, and the structures grow according to the linear theory. On small scales, 
the gravitational collapse is non-linear, and the accurate behaviour of the density contrast 
can only be computed by using N-body simulations. 

Besides its intrinsic phenomenological character related to the large scale structures, the 
correct interpretation of the matter power spectrum is concerned with one of the experimen- 
tal techniques for determining the mass of neutrinos through cosmological measurements, 
namely the CMB results. It is done through the inference of a transfer function in the 
matter power spectrum at small scales [7, 8] . The contribution due to massive neutrinos to 
the closure fraction of cold dark matter at present can be depicted from the modifications 
on the matter power spectrum, even for neutrinos behaving like hot dark matter at higher 
redshifts. It follows that the amount of cold dark matter at earlier epochs should be substan- 
tially reduced, suppressing the formation of large scale structures, when it is compared to a 
situation without massive neutrinos [9-11]. The effects, however, are attenuated in case of 
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hot dark matter or massless neutrinos. A comparison of the observationally-inferred matter 
power spectrum with the power spectrum expected without the effects of massive neutrinos 
allows one to compute the effective contribution of neutrinos to the cosmic inventory. 

The point at our manuscript is that the hypothesis of a tiny fraction of the cosmic 
inventory evolving cosmologically as a degenerate Fermi gas (DFG) test-fluid can mimic 
the massive neutrino cosmological behaviour at several dominant cosmological backgrounds. 
The approach for treating the cosmological perturbation for ultra-relativistic neutrinos as 
hot dark matter, in general, differs from that for treating non-relativistic neutrinos as an 
additional cold dark matter component. Treating massive neutrinos as a DFG not only gives 
novel ingredients to the recursive idea of a self-gravitating Fermi gas model but also allows for 
quantifying the smooth transition between ultra-relativistic (UR) and non-relativistic (NR) 
thermodynamic regimes for neutrinos in the cosmic background. In a previous issue [12], an 
analytical procedure based on Bianchi identities was adopted for obtaining the evolution of 
perturbations for a class of fluids which evolve from relativistic to non-relativistic regimes, 
a supposition that is usually adopted for cosmological neutrinos. Our proposal becomes a 
convenient tool for describing the role of such particles in the cosmic inventory so that some 
peculiar modifications on the matter power spectrum can be debugged. 

Assuming that the neutrinos temperature at present should be of about Tq^, ~ 10~^K 
(depending on some phenomenological correspondence with the neutrino mass), and com- 
paring this with a tiny fraction of a DFG test-fluid [12] in equilibrium with the radiation 
background in the beginning of the radiation-dominated era [13], one should expect that the 
DFG's contribution to the matter power spectrum could be of the same order of magnitude 
of the neutrino's contribution. 

Besides theoretical speculations about the existence of some kind of DFG produced in 
the early universe [14-16], and the idea of a self-gravitating Fermi gas model introduced to 
explain the puzzling nature of white dwarf stars [17], it may be expected that at some stage of 
the evolution of the Universe, primordial density fluctuations have become gravitationally 
unstable forming dense lumps of dark matter. Such possibilities of forming large scale 
structures thus stimulate our meticulous investigation of the behaviour of a DFG as a test- 
fiuid in the cosmological background of radiation, matter and cosmological constant (A). 
Despite not necessarily resulting in matter lumps, we shall notice that it somehow contributes 
to the matter power spectrum related to large scale structures. 
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To quantify some analytical aspects related to the here proposed DFG approach, which 
results in fiducially modified and numerically confirmed predictions for the matter power 
spectrum of the cosmic inventory, our manuscript was organized as follows. In section II 
we report about the formalism of linear perturbations focused on the study of a perfect 
fluid through the Einstein equation. In section III we obtain some analytical expressions 
for the density fluctuation, S, and for the fluid velocity divergence, 9, in case of a DFG 
at radiation-dominated, matter-dominated, and A-dominated background scenarios, and for 
super horizon solutions in case of a radiation-to-matter transitory regime. In particular, we 
consider the conservation equation for the stress-energy tensor with a vanishing anisotropic 
stress tensor (a). Reporting about the systematic way of accounting for all the cosmological 
perturbations, i. e. through the Boltzmann equation, in section IV we quantify the be- 
haviour of the inhomogeneities in the Fermi gas at the radiation-dominated era, through an 
ultra-relativistic analytical approach, and at the matter-dominated and A-dominated eras, 
through a non-relativistic analytical approach. All the analytical results are consistently 
verified through the numerical calculations described in section V, where the DFG contribu- 
tion to the matter power spectrum is compared with previously obtained numerical solutions 
for massive neutrinos corresponding to a tiny fraction of the dark matter components. We 
draw our conclusions in section VI. 

II. COSMOLOGICAL BACKGROUND AND LINEAR PERTURBATIONS 

The cosmological evolution of a homogeneous Friedmann-Robertson- Walker (FRW) fiat 
universe with averaged energy density, p{ri), and averaged pressure, P{r]), is described in 
terms of the scale factor, 0(77), through the following components of the Einstein equation. 



where 77 is the conformal time, G is the Newtonian constant, and we have set the light 
velocity equals to unity, i. e. c — 1. 

From the properties of a perfect fiuid, the above elementary textbook equations can be 
depicted from the conservation properties of a stress-energy tensor that, at a comoving 
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frame, is given by [10] 

rj^i^i^ ^ Pgnu + (p + r) U^'W , (3) 

where U is the 4-velocity of the fluid particles, g^i, is the metric tensor, and energy den- 
sity and pressure are decomposed into averaged and perturbative values as p — p + Sp and 
V — P + SV. Once one assumes the vahdity of the Einstein equation, the corresponding con- 
servation equation for a non-interacting single-particle fluid should be given in terms of the 
covariant derivative, T.^" = 0. One shall notice that the collective behaviour prescribed by 
our DFG approach guarantees that the stress-tensor is isotropic and diagonal, and that the 
energy-momentum tensor can thus be analytically described in terms of the abovementioned 
space-time dependence. 

Throughout our study, we have chosen to work out the perturbation equations for the cos- 
mological scenario in the longitudinal gauge since it has set simpler and clearer connections 
to preliminary studies on cosmological perturbations [7, 9-11, 18]. By considering that the 
non-diagonal metric perturbations vanish [11], two scalar degrees of freedom are eliminated 
from Eq. (3), and the remaining ones are defined in terms of two scalar potentials, 4> and ip, 
as 

g^o^-a\l + 2il;)5^o, ft,- = a'(l - 20)5^ . (4) 

The Bianchi identities [10] then provide us with the constraints on the stress-energy tensor 
in the Fourier /c-space, which can be summarized by the continuity equation, 

5^-{l+uj){e-3^)-3-{l + 5){l + uj) - ^{1 + 5), (5) 

a p 

and by the Euler equation, 

9 = ; ak"^ - A-9 + ibk'' - , (6) 

l + u 5p a ^ \ p l+u' ^ ^ 



with 



dp 

a— + 3p{l + uj) = 0, (7) 



where uj — V/p, and 9, S and a are, respectively, the density fluctuation, the fluid velocity 
divergence, and the shear stress deflned by 

S = Sp/p, {p + P)9 = ik^ST^, and {p + P)a = - (^k% - ^S'^ E^, (8) 
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with Ti^Ti- {5ll?>)Tl 

The solutions of Eq. (1) for radiation- and matter-dominated universes, and for the 
cosmological constant, with the corresponding equation of state respectively represented by 
V-y = P'^l'i-, Vm = and Pk = —pA: are given by 



- 

Pm — PO — ^ 



and Pa ^ pqQa, (9) 
where the scale factor during an approach for radiation-to-matter transition era is given by. 



STiGpo 



(10) 



and during the A-dominated era by 

)l/3 



a(r]) 



1/3 



SirGpoil/ 



(11) 



where po is the Universe's density at present, and = Ps/ Po, with s — ^, m, A. 

Finally, by decomposing the Einstein equation into time-time, longitudinal time-space, 
trace space-space, and longitudinal traceless space-space components, one obtains four equa- 
tions to the linear perturbations in the /c-space, which through Eq. (4), in terms of and 
ip, are given by 
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(12a) 
(12b) 

(12c) 
(12d) 



l2T:Ga'^{p + V)(T^ A;2(0- V). 

In the limit of "single-particle" fluid eras separately described by each element of Eq. (9), 
and in the absence of the shear stress, i. e. with o" ~ 0, one can set -0 ~ in order to 
derive some analytical expressions for the evolution of the metric perturbation through scale 
factor dependence on 77 given by Eqs. (10) and (11). The full continuity, Euler and Einstein 
equations thus result in 



, / \ ^ cos uj + uj sm uj ^ sm uj — uj cos uj 
(t^iVn) = Ci VC2- 



(13a) 



Mv) = -^ + D2, (13b) 

Mv)-Hb-Vf + F2{b-V), (13c) 

with u! = k?]/ \/3. The coefficients with sub-index 1 in the above equations are related to the 
so-called iscntropic decaying modes, in the same fashion that the coefficients with sub-index 
2 are related to the isentropic growing modes. As expected, decaying modes are suppressed 
during the cosmological evolution. In addition, the solution for at the super-horizon limit 
during the radiation-to-matter transitory period is described by 

'l>iy) = Si—^ + S2 , (14) 

where y = a/aeq = Pm/p^,, a^q is the scale factor at the time of equality, i. e. when 
and a = 0(77) is given by Eq. (10). 

III. DFG TEST-FLUID SOLUTIONS 

To obtain analytical solutions for a DFG as a test-fluid in each one of the cosmological 
background scenarios briefly discussed in the previous section, one has to follow a sequence of 
consistent simplifications. Let us consider that the phase space distribution of the particles 
gives the number of particles in a difTerential volume dx^dx^dx^dpidp2dp^ in phase space, 

f{x\ pj, ri) dx^dx^dx^dpidp2dps/{2nf = dN, (15) 

with Pj = Pjip, rij), and where / is a Lorentz scalar that is invariant under canonical trans- 
formations. The general expression for the stress-energy tensor, in case of a single-particle 
fluid, written in terms of the distribution function / and of the physical quadrimomentum 
components, P^, is then given by 

^'^ " j (27r)3(-^)i/2 po /(^ ^PJ^^I) ' (16) 
The phase-space distribution function evolves according to the Boltzmann equation as 

Df df ^ dx' df ^dqdf ^ dn,, df / df\ 
drj drj drj dx^ drj dq drj dui \ drj j ' 

where, for non-interacting fluids, the right-hand (collision) term vanishes. The fluid 
temperature is, however, not uniform. Thus one conveniently sets T — T + 5T with 



ST — 5T{x\nj,r]), ap^ = = qn\ where q is the comoving momentum, and the time- 
dependence is relegated solely to T. Expanding the distribution function, f{x\pj(p,nj),'r]) 
up to first-order in the temperature perturbation, ST, with f\f = f{p,T{r])), one has 



5/ 
dT 



6T{x\nj,7]) = f{p,7]) [1 + ^{x\pj,7])] , 



(18) 



and the Boltzmann equation in the Fourier /c-space for the longitudinal gauge becomes 



(1 + + + i^{k ■ n)^f + -ie{k- n)i^j 



drj 



dq 




(19) 



where energy and momentum are respectively rewritten in terms of the (pseudo)comoving 
energy, e = aE = a{p^ -\- m'^Y^'^ = (g^ -|- a^m^Y^"^ , and of the comoving momentum, q = pa. 
By reducing the notation for the averaged temperature, T, to ~ T, the distribution function 
for a DFG can thus be given by 



f{p,v) = (exp 



E- fi 
T 



+ 1 



gd for E < II 
ioi E > II 



(20) 



where ga is the number of spin degrees of freedom, ix is the chemical potential (that for 

lower temperatures approximates the Fermi energy ix > Ep), and Planck and Boltzmann 

constants, h and /c^, were set equal to unity. Reporting about Eq. (16), the pressure and 
the energy density for a DFG are given by 



=9d 



Pd =9d 



m 



167r2 



m 



-X 



In 



1 + 



X 



x" 



(2 + x')-ln 



1 + 



X 



(21) 



where % = {m/qp)a, m is the particle mass, and qp is the comoving Fermi momentum. 

Assuming the fluid behaves as a test-fluid under perturbations due to a background scalar 
potential, 0, the corresponding continuity and Euler equations would be given by 



d5pa 
df] 



m 



37r2 x' 



30 



« 4 + 3x' 
dpd- 



a 



1 + x' ' 



a 

Od + crdk'^ + -— — ^ -(j)k'^ ^ 5pa 
al -Fx 



m4(l + x2)3/2 



(22a) 



(22b) 



where we have identified the DFG through the index d. Supposing that the DFG averaged 
density is equivalent to the averaged density of massive neutrinos, and considering that 
neutrinos are ultra-relativistic only deep inside the radiation-dominated era, one can simplify 

the above obtained equations by setting x ^ 1 during the radiation-dominated era (when 
SVd/Spd ~ Vd/ pd ~ 1/3) and by x > 1 after the radiation-to-matter transitory regime (when 
^Vd/^pd ~ l/3x^ and Pd/pd ~ l/5x^), which is quantified by a = Oeg. 



A. Radiation-to-matter transitory regime at very large scales 



For krj and with vanishing u, Eq. (22) becomes 

din [x{pd + Vd)ed] 



dhax 



9^ m^Jl + ^9(l) 

X—opd - -j^ — — ^ + ^Pd 



= -3, 
4 + 3x' 



<9x TT^ x" dx 1 + 



0. 



(23a) 



(23b) 



Using the analytical solution for in case of super- horizon scales (c. f. Eq. (14)), one thus 
obtains 



71 



m 



5d = ^ 



y/TT^ 16 + 8y-2y^ ' 

Ol 5 r 02 n r <J4 



(2 + X') 



larccschx 



(24) 



(25) 



yO 

where Si, with i = 0, 1, 2, 3, are constants used to independently match preliminarily un- 
coupled solutions for radiation- and matter-dominated eras. Although the density contrasts 
for cold dark matter (CDM) and photons remain the same for the period of radiation-to- 
matter transition, Fig. 1 shows that the density contrast for a DFG does not reach the same 
growing rate of that obtained for CDM at the end of the matter-dominated era, in case 
of super-horizon scales. We shall see that such a behaviour recurrently happens at scales 



> 



Mpc/h, with a non-trivial contributions to the matter power spectrum today. 



B. Radiation-dominated era 

For the radiation-dominated era with vanishing cr, the components of Eq. (22) (or Eq. (5) 
and Eq. (6)) can be approximated by 

5d^-^ed + ^, and Od^ k^^5d+ k''(l>. (26) 



Using the analytical solution for from Eq. (13a), one thus obtains 

2(1 — uj'^)cosuj + u!(2 — uj'^)smu} 2(1 — o^^) sino^ — a;(2 — a;^) cosa; 
Od = 2Ci \-2C2 ^ , (27a, 




(27b) 



where Cj, with i = 1, 2, are also fitting constants. In this case, the DFG perturbation 
reproduces the same behaviour of radiation perturbations since, for the ultra-relativistic 
limit, one considers x ^ 1- The solutions for the growing mode related to the coefficient C2 
correspond to oscillating perturbations. It occurs because the non-vanishing pressure has 
prejudiced the gravitational collapse. 



C. Matter-dominated era 



Deep inside the matter-dominated era, when the critical density approximates p^, also 
with vanishing a, the components of Eq. (22) become 

7 2 X 

Sd^-9d + 3<i>, and 9a ^ ^ - k'' aa+ ^ ^ - % . 
Using the analytical solution for from Eq. (13b), one finally obtains 

7.2 „2 



(28) 



12V3c; 



^-12^ + 96c: 
5 r)^ 



-A 



and 



V^k ^ V^k VSk ^ VSk 
cos Ci- h sm Si- 



k^T]' „ VSk ^ . Vsk 
h Us COS U4 sm 

6 2Cm Tj 2Cm Tj 



(29) 



^-^ + ^ + 256c; 
Tj'^ 15 



+D2 



8c3 



"m 



Ci 



77^ 

. V^k . y^k ^. V^k ^/3k 
1 sm bi cos 

1] T] 2Cfn Tj 2Cm f] 



V^k . V^k ^ s/2>k V^k 

sm Da cos 



2cmr) 



k^ 2r]' 



,2 ^2 



2cmr) 

1 

(30) 



where = {G pQrm:VL„i) / qp^ and Si(2;) and C\{z) arc the integral sine and cosine func- 
tions. The term proportional to at Eq. (29) shows that the density contrast, 5, grows 
at the same rate as the expansion parameter, a. The last term between brackets oscillates as 
sin(-\/3A;/(2cm ?7)) for 77 < y/3k/{2cm), after which, it assumes a logarithmic growing depen- 
dence on the conformal time, 77. Turning to the solution for the fluid velocity divergence, 6, 
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described by Eq. (30), the terms inside the first brackets describe damped oscillations which 
are similar to those parameterized by sm{^/3k/{2cJnrj)) at Eq. (29). For rj < y/?>k/{2cm), 
the decreasing amplitude is proportional to the inverse of the expansion parameter, i. e. to 
1/a, after which the oscillating behaviour is completely suppressed. The last term in the 
brackets presents a linear growing dependence on the conformal time, i. e. 77 oc ^/a. 

In general lines, the perturbations described by Eqs. (29) and (30) reproduce an oscillating 
behaviour at the beginning of the matter-dominated era up to the point where it can be 
suppressed by the subsequent growing modes. The complete solutions for 5 and 9 can be 
depicted from Figs. 2 and 3, from which one can notice that DFG perturbations never reach 
the same growing rate of CDM perturbations. 



D. A-dominated era 



If one assumes a non-relativistic dynamics for the DFG test-fluid at late times, that 
happens when one sets x » 1 at Eq. (21), the cosmological constant era can be reproduced 
by the same dynamics that drives the perturbations into the matter-dominated era, (c. f. 
Eq. (28) with a given by Eq. (10)). Using the analytical solution for (/)a from Eq. (13c), one 
thus obtains the density contrast 



5d=3Fi4(6-77) + F3Cos 



TT 77(26 — rj) 
2 d 



+ 3ca{F2 + Fi4}|cos [| 



+ 



TT 



cos 



z 

TT 



-|- F4 sin 
Fc [z\ -\- sin 
Fs [z] — sin 



TT 



TT ri{2h — rj) 
2-^— 

Fs[^]} 
FcN} , 



(31) 



7r 



where z = (b — r])/c\, c\ = \/3cm7r/kqF, = cm/qp, Fs{z) and Fc(^) are respectively 
the Fresnel integral sine and cosine functions, and c and b are constants respectively given 

by 

1 



4a 



l/2> 



(32) 
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In correspondence with the above density contrast, the solution for the fluid velocity diver- 
gence is given by 



n 



9 Fi{ri - by + F^i—z cos 
Ca 



TT r]{2b — rj) 
2 d 



cos 



+ < COS 



12 . 


Fc [z] 


+ sin 




Fs[z]] 




Fs [z] 


— sin 




Fc[z]] 



P4 — 2;sin 
Ca 



3tt 



TT 77(26 — 77) 

2^r~ 



(33) 



F2 — + Fi"\/3/ccaC;^ \ ct^z 
Ca 



These solutions correspond to decreasing functions that result on decaying modes of pertur- 
bation, an effect that match the dynamics ruled by the negative pressure of the dominant 
cosmological background. 



IV. BOLTZMANN EQUATION SOLUTIONS FOR A DFG 

In certain sense, the thermodynamics of a DFG is driven by the cosmological behaviour 
of the chemical potential, ~ /i(a). However, depending on its effective contribution to the 
fluid dynamics, which sometimes is not relevant, additional contributions to the analytical 
characteristic of the distribution function are usually discarded by setting ^ 0. That is 
not the case of our approach, since we are interested in depicting the cosmological behaviour 
of a DFG. The simplest way to consistently include the chemical potential in the following 
calculations is through the zeroth-order Boltzmann equation for a non-interacting fluid (c. 
f. Eq. (19)) as 

from which the last equality sets the constraint among the energy, E, the temperature, T, 
and the chemical potential, fi. 

The condition given by the above equation results in the following constraint for the 
dependence of the distribution function on the temperature. 



9/ 
dT 



_ aE n- Edf _ e u -edf 

T 

with p = a/i, a kind of comoving parametrization of the chemical potential. It allows one to 
identify the first-order perturbation coefficient, ^, that appears into Eq. (18), as 

^ix\nj,r)) = A_l — A{x\nj,r)), (36) 

12 



where A = 5T/T. By substituting the above expression for ^ into Eq. (19), the Boltzmann 
equation can be reduced to 

where A{x\nj,r)) is read from V{k\nj,rj) in the Fourier /c-space as 

A{x\ rij, T]) = J dk^ exp [ikix'] V{k\ Uj, r]). (38) 

The dimensionahty of the problem can be reduced by noticing that the evolution of Eq.(37) 
depends on the direction fi parameterized by the angle related to A; • n = cos 93, which is a 
natural consequence of the isotropy of the homogeneous background. 

Performing the P/-Legendre expansion with respect to cos((/?), one obtains 

00 

V{k\n^, ri) = ^(-i)'(2Z + 1) A^^, v) Pi{cos{^)) , (39) 
1=0 

that results in 

Ao = ^ J dflV, Ai = ^ J dfl cos^V and A2 = --^ J dfl (^cos^ ^ - ^ V , (40) 

for the first three multipole cofficients. The corresponding evolution equations for the mul- 
tipole coefficients, Ai{k,ri), are obtained through the integration of the Boltzmann equation 
(37) multiplied by Pi{k ■ h) = Pi{cos (f) over the solid angle dfl = d9 d{cos if), that results in 

d le{u-e)dl \ ^ df df 
or] \ q oq j oq oq 

d f^(„_e)df \ df k 

^(^^^Aj+(.-e)-— [(/ + 1)A.,,-/A._,]^0, for l>2. 

(41) 
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The perturbations for the stress-energy tensor from Eq. (16) can thus be rewritten as 

'^^'^^ ^^^^ ("-^^a^' 



67r2a4 J dq ' ^^2) 

Od{Pd + Vd)=^^^ J dqq\{u-e) — , 

adiPd + Vd)--^, jdqq i^-e)^^, 

and finally, after integrating Eqs. (41) over q and using the results from Eq. (42), one 
can resolve the system of coupled equations and find the time dependence of the DFG 
perturbations. From such results we shall verify that it is possible to analytically depict the 
DFG transitory regimes, from relativistic to non-relativistic dynamics. 



A. Ultra-relativistic DFG during the radiation-dominated era 

When particles are ultra-relativistic, i. e. e — q, and u is set equal to zero, the above 
results should be reduced to that obtained for massless particles [10, 11]. Differently from 
which is observed when one sets the ultra-relativistic approximation to /, the spectrum 
does not remain Planckian when 1/ ^ 0. The g-integration can be performed separately 
and the dependence on x\ rij and r] is kept active. The cosmological behaviour of the 
temperature depicted from Eq. (17) is given by T{rj) ~ Tom/[a{ri) qp] since the energy, E, 
is approximated by its ultra-relativistic dependence on the momentum (~ p) through the 
relation 

— 7^ — ~ Qf , 43 

where Tq is the fluid temperature at present. By performing the quoted integrations for 
Eqs. (41), the multipole evolution equations become 

Ao + kAi-^cp^O, (44) 
Ai-^(Ao-2A2)-^0Ci. = O, (45) 

Ai + ^^[{l + l)Ai+i-lAi_i]^0, for l>2, (46) 
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and 



where cp = q%/{mTQ\n{2)). By considering the analytical solution for (f)^ from Eq. (13a), 
and assuming that the contribution from higher order multipoles A/, with I > 2, are sup- 
pressed, the above equations lead to 

C c C c 
Aq — — 5— [cos a; (l — 00"^) -|-a;sina;lH — [sin a; (l — a;^) — a;cosa;l-|-C3COsa;— C4sina; , 

(47) 

and 

CiCp , . 1 C2CF p . 1 , C's . , C4 /.qx 

Ai = p — [cos oj + UJ sm uj\ p — [sm u — u cos u!\ H — ^ sm a; H — p cos a; . (48j 

^/3u!'^ v3a;2 V3 v3 

The above results can be substituted into the integrals from Eqs. (42) in order to give 

5,«2Vln(2)^Ao, (49) 

~ ^ (501 

rri 

^df:al27rMn(2)^A;Ai, (51) 

and, eventually, 

a, ^ -2V ln(2) ^ A^, (52) 

that reproduce the results from Eqs. (27a) and (27b) obtained through the continuity and 
Euler equations, i. e. Eqs. (29) and (30), once one has assumed that the DFG behaves as 
an ideal fluid. 



B. Non relativistic DFG during the matter-dominated era 

Deep inside the matter-dominated era, when ~ Pm, the cosmological behaviour of 
the temperature depicted from Eq. (17) is given by T{r]) ~ TQa{r])~'^ since the energy, E, 
is approximated by its non- relativistic dependence on the momentum (~ /{2m)) through 
the relation 

By evaluating the quoted integrations of Eqs. (41) over the momentum 5, and assuming that 
A; ~ 0, for Z > 2, the multipole evolution equations can be written as 

Ao + - Ai - cf4> = 0, (54) 
X 
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and 



k 



A,-^ (Ao-2A2)-^XCF0 = O, 
3x 3 



[(Z + l)A^+i-ZA,_i]=0, 



for I > 2. 



(55) 



(56) 



.2/ + i; X 

By substituting the analytical solution for 0^ from Eq. (13b) into Eqs. (54-55), one obtains 



Ao 



D2 Cf 



k^ 



— - 12— + 96c^ 

5 rf ^ 



cos Li h sm 



\flk ^ . VSk 
+ I?3 COS D4 sm 



2cm ?7 2c^ 7] 



-Si- 



2c^ 2cjn Tj 



D2 Cf 



2C}72 ?7 (5T) 
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Ai =DiCf 
+D2 Cf 



y/SCrn 

kyvs 



96A;2c2 
rf 



+ 



6A;4 



^. VSA; . \f2>k 
Ci- sm 



15 



+ 256c^ 



sm ^ cos 



^. ^A; V3A; 

bi- cos- 

2cm-n 2c^ 77 2c^?7 2c„?7 



V3 2c^77 V3 

kr] '"^ 



+ D2CF 



9c 



y3A; 
2c,n77 

^ , 2^ 
4 3 



(58) 



where Cm = {G pomn flm)/(lF, as it was previously defined. Evaluating the integrations 
prescribed by Eqs. (42) and using the non-relativistic approximation for the DFG through 
Eq. (53), the perturbations for the stress-energy tensor are given by 



127r3ln(2)^An 



1 



5pd 3x 



2' 



Od ~ 127rMn(2) 



mTo Ai 



and, when required. 



<7,«87rMn(2)^^. 

Tf X 



(59) 

(60) 
(61) 

(62) 



Again, the above solutions reproduce the same behaviour of the solutions from Eqs. (29) 
and (30) for an ideal fluid approach. One should notice that the non-null perturbation on 

the averaged pressure, as obtained through the ideal fluid approach, explains some small 
deviations from pure cold matter perturbations, for which the pressure should vanish. 



C. Non-relativistic DFG during the A-dominated era 



As we have stated in the previous section, when one assumes the non-relativistic dynam- 
ics for a DFG test-fluid at late times by setting % » 1 at Eq. (21), the dynamics for the 

16 



cosmological constant era can be reproduced by the same equations of matter perturba- 
tions. In this case, it is reproduced by Eqs. (54-56). By following the same procedure as 
for radiation- and matter-dominated eras, i. e. using the analytical solution for 0a from 
Eq. (13c) into Eqs. (54-56), with A; ~ 0, for / > 2, one obtains 

TT r]{2b — T]) 



Ao ^F^Cf 4(1] -b) + Fs cos 



TT 77(26 — 7]) 

2 d 



+ Fi sin 



+ 


|cos 




c [z] + sin 


"TT 2" 

—z 
12 J 


+ 


l^cos 




s [z] — sin 


[r1 



cf 



kir „ Q 3 

71 -^"^^ 



Ai = - Fi Cf cl — {b -rj) + ^ cos 
TT V3 



TT 77(26 — 77) 

2 d 



sm 



TT 77(26 — 77) 
2 d 



— |cos 


2I 

.r. 


Fc [z] + sin 


.r. 


-|- |cos 


\^ 21 

.r. 


Fs [z] — sin 


\^ 21 
.r. 



Fs \ z]\ If, ci/' 



2^X ^-^l4- 



Fc[z]}{F, + F,cl} 



CpCA 



(63) 



(64) 



with c\ = y/3c'mn/{kqF), = cm/qF, where c and b are given by Eq. (32). 

Since the Eqs. (59-62) are still vahd, the DFG perturbations during the A-dominated 
era (obtained using the Boltzmann equation) are equivalent to the solutions obtained using 
Eqs. (31) and (33). 



V. NUMERICAL SOLUTION 



Herewith we shall consider the numerical procedures for treating perturbations valid for 
test-fluids in a DFG regime. For such a purpose we have compiled the pubhc code CAME 
[19] based on the synchronous gauge scalar perturbations on the flat Friedmann-Robertson- 
Walker metric with isentropic initial conditions for the stress-energy perturbations. As one 
shall notice, the numerical results that we have obtained are in complete agreement with 
the preliminary analytical studies that we have proposed in the previous sections. 

Considering that the perturbed flat Friedmann-Robertson- Walker metric in the syn- 
chronous gauge is described by 

ds^ = [-^77^ + {Sij + hij) dx'dx^] , (65) 
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the continuity and Euler equations of the stress-energy tensor for ideal fluids can be written 
as 

d5p I /) a\ ci f 

^=-p-(p + P)|^«+- + 3-j-3^.p(l + ^), (66) 

and 

e + ak'' +4^^^ {p + P)+e(p + P^ ^5pk^^-^, (67) 
with h = ha. In analogy to Eq. (21), one thus obtains 

and 

which are used in the numerical calculations We have performed the numerical calculations 
to obtain the time dependence of the DFG perturbation by introducing the constraints given 
by Eq.(68) into the CAMB code, with 0. 



A. Background and perturbative initial conditions 

To compare DFG with massive neutrino perturbations in the background of the ACDM 
model, we have set coincident values for the averaged densities of both fluids as preliminary 
conditions, at early times (radiation-dominated) and at very late times (matter-dominated 
era). To obtain the comoving Fermi momentum qp for a DFG, we have considered that the 
ultra- relativistic hmit of Eq. (21), 

MX«1)«3^^ = |^|, (69) 

does not depend on the mass, m. By equating the averaged densities of the equivalent matter 
components: DFG and neutrinos, i. c. by setting pd = (assuming three degenerated 
families for each one of them) , one easily verifies that 

7/4x4/3 

Pd{x < 1) = Pu{a < Oeg) = 3g ( 1 P^{a) , (70) 
which results in qp ~ 3.65 • 10~^ eV/c. 

18 



By introducing the obtained value into the non-relativistic approximations for the DFG 
averaged density, 

Mx»l)-^^ = ^^-, (71) 
and following the same equality, pd — Pv, in the limit where a tends to unity at present, 

Pdix > 1) = ^ 1) = Po% , (72) 

one can easily deduce the value for an effective neutrino mass given by m ~ 0.24 eV/c^. 

Finally, by integrating the isentropic initial conditions over q for the multipole expansion 
of the massive neutrino perturbations [10], and using the DFG phase-space distribution, we 
reproduce the same initial conditions for DFG perturbations as those for massless neutrinos. 

Fig. 4 shows the density contrast, S, that we have obtained for a DFG. During the 
radiation-dominated era (a <IC Og ~ 2,7 • 10~^) the DFG density contrast, Sd, has the same 
oscillating behaviour of radiation perturbations, 5-y, since both fluids have approximately the 
same evolution equations. It occurs for a DFG approximated by ultra-relativistic conditions. 
During the matter-dominated era the DFG density contrast oscillates up to the beginning 
of the non-relativistic regime. After being transposed to the non-relativistic scenario, the 
DFG density contrast oscillating modes are suppressed and the growing modes reproduce 
the same behaviour of the CDM perturbation. Finally, when the Universe passes to the 
A-dominated era, the DFG density contrast growing modes are suppressed by the negative 
pressure of the dominant fluid. 

Fig. 5 allows us to depict the difference between the density contrast of a DFG, for 
which we have suppressed the higher order multipoles, and the density contrast of massive 
neutrinos evaluated numerically. The time dependence of the perturbations for both fluids 
has two fundamental differences. Firstly, as a vanishes only for our DFG model, the damped 
oscillations for the massive neutrino density contrast are not observed for the DFG during 
the radiation-dominated era (ultra-relativistic period). In addition, since we have followed 
the above approach for matching DFG and neutrino energy densities, the neutrino mass 
for a completely non-relativistic fluid is supposed to be given by m,^ ~ 0.46 eV, while 
the DFG was already stated as rrid ~ 0.24 eV. Due to such a peculiarity, the massive 
neutrinos become non-relativistic earlier. These two opposed effects compensated themselves 
so that the density contrast for massive neutrinos and for a DFG have equivalent amplitudes 
at present. The matter power spectrum for hCDM -\- v and KCDM -\- DFG universes 
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are therefore very similar for large and small scales, with a small relative deviation for 
intermediate scales. 

Fig. 6 shows the difference between the matter power spectrum, P^(fc), for the cosmic 
inventories described at Table I. The matter power spectrum is defined in the Fourier /c-space 
as the averaged of the correlation function of the density contrast as 

{5l{k, 7l)6m{k', 71)) = {27r)'Pm{k) S'{k - k'). (73) 

By observing that 5m = Pi^i) / Pi) (^^^ ^ matter components) and integrating the 
above equation over the solid angle dQk, one can depict an explicit value for Pm{k). All 
the resulting matter power spectra were normalized to match one each other at large scales. 
Fig. 6 shows a certain level of decreasing on the ACDM power spectrum at small scales. This 
is due to the introduction of an additional component into the cosmic inventory, namely 
massless neutrinos, massive neutrinos or even a DFG test-fluid. Such a large difference 
is indeed evident for small scales, i. e. when one compares the density contrast between 
CDM plus barions and CDM plus barions modified by additional (and eventually, exotic) 
components, as it has been discussed in some previous references [11]. 
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TABLE I: Cosmic inventory corresponding to the matter power spectrum of Fig. 6. It is assumed 
a flat universe withSjOj = 1 that results in r = 0.027454 

Just to end up, in what concerns the neutrino degeneracy, one can ascertain that the 
value of the chemical potential of electron neutrinos (/Xe/To) should be close to —1 in order 
to explain with an acceptable confidence level the phenomenology related to the observed 
variation of deuterium by roughly an order of magnitude [8]. It would induce a variation 
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in the total energy density during the radiation-dominated stage, which is excluded by the 
smoothness of the CMB. However, this objection could be avoided if there was a coincidence 
between different leptonic chemical potentials such that in different spatial regions they 
could have the same values but with the interchange of electronic, muonic and/or tauonic 
chemical potentials [20]. For further details and discussions on the energy dependence and 
effective chemical potential one might address to preliminary calculations [20, 21] where the 
main points are that they are not obtrusive to our results. 



VI. CONCLUSIONS 



The hypothesis of a small fraction of matter components evolving cosmologically as a DFG 
test-fluid in some dominant cosmological background was evaluated through numerical and 
analytical approaches. Our proposal was based on the fact that for any type of fermionic dark 
matter, which includes relativistic and non-relativistic neutrinos, the average phase-space 
density cannot exceed the phase-space density of the DFG, which leads to the possibihty 
of forming extended overdense regions of a fluid of degenerate particles. In the scope of 
our analysis, we have introduced some analytical procedures for obtaining the evolution of 
perturbations for some categories of relativistic and non-relativistic test-fiuids. Expressions 
for the time evolution of density contrasts and fluid velocity divergencies with vanishing 
anisotropic stress, o", were obtained for a gas of massive and degenerate fermionic particles 
in the background regimes of radiation-, matter- and A-dominated universes, and for super- 
horizon scales during the radiation-to-matter transition period. Mathematically, one can 
depict from the set of Eqs. (42) that the anisotropic stress is proportional to the quadrupole 
term in the same fashion that the DFG pressure perturbation is proportional to the monopole 
term. In this case, the same arguments that one could use for disregarding the contribution 
from cr to the evolution of the perturbations for ultra-relativistic and non-relativistic are 
analogously maintained. In case of cosmic neutrinos, one certainly should not discard the 
contribution due to the anisotropic stress (multipole Z = 2) in case of analyzing GMB and 
CuB anisotropics due to scalar perturbations and the Gi^B coupling to gravitational waves 
due to tensorial perturbations [22, 23]. Evidently, that was not the case of the analysis 
that we have performed. Anyway, in case of cosmological massive neutrinos, the collision 
term is suppressed since the collision rate succumbs under the Hubble rate, and neutrinos 
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free-stream. Specifically concerning the contribution to the matter power spectrum, which 
comes significatively from non-relativistic neutrinos (during matter-dominated and matter- 
to-A present epochs) the I > 2 multipole contribution are completely suppressed. 

Turning back to our primary results, during the period when the DFG could be ap- 
proximated by an ultra-relativistic gas, it was obtained that the density contrast presents 
an oscillating behaviour similar to that of the photon perturbation during the radiation- 
dominated era. Through analogous comparison with matter components, the freezing of 
DFG components leads to a combined linear plus logarithmic growing mode for the density 
contrast, similar to which one could obtain for CDM perturbations. Finally, during the 
A-dominatcd era, the DFG perturbation terms are suppressed (as 1/a) due to the negative 
pressure of the dominant fluid. 

Through numerical calculations, we have quantified the impact of a DFG component into 
the cosmic inventory on predictions about the matter power spectrum. It was performed by 
running the public code GAMB based on the synchronous gauge and comparing the results 
with the previous analytical studies based on the longitudinal gauge. For intermediate 
regimes, one could observe the suppression of the growing adiabatic modes which turns 
into a smoothly oscillatory behaviour. Also for the fluid velocity divergence, the transition 
from ultra-relativistic to non-relativistic regimes was quantified. In this sense, our results 
have become a convenient tool for performing preliminary tests for HDM to GDM test-fiuid 
configurations. Being more specific, the effects quantified through analytical and numerical 
calculations and their relative impact on the cosmological structure formation were obtained 
by substituting the massive neutrino fiuid by a DFG matter component. Assuming isentropic 
initial conditions, the Fermi momentum, qp, and the mass value, m, were adjusted in order 
to give early (radiation-dominated) and late (matter-dominated) times averaged densities 
that reproduce the cosmological background neutrino densities obtained through standard 
predictions from the AC DM model. 

Generically speaking, the global evolution of the DFG perturbations obtained by numeri- 
cal procedures is consistent with the results interpreted from analytical solutions. During the 
ultra-relativistic regime, the massive neutrino density contrast presents a damped oscillation 
while the power is transferred to higher multipoles. Obviously, through the approximation 
used for massive neutrinos, they become non-relativistic earlier than a DFG since m < rrii, 
for the same averaged densities. Meanwhile, the density contrast for the DFG and for mas- 
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sive neutrinos are very similar during the non-relativistic regime (c. f. Fig. 5), even though 
it does not reach the growing rate of the CDM density contrast. The effects on the power 
spectrum for large and small scales (A;?7 <^ 1 and krj 1 ) are minimal if changing massive 
neutrinos by DFG. One of the subtle points of our analysis is that the effects on the matter 
power spectrum, due to the conversion of massive neutrinos into a DFG fluid in the cosmic 
inventory, are relevant only for intermediary scales as depicted from Fig. 6. 

To summarize, the reproducibility of our results through the conservation equations de- 
rived from the Bianchi identities, through the analytical approximations on the Boltzmann 
equation, and through the numerical calculations, ratifies the consistency of our analysis. 
Once the region of intermediate scales for the matter power spectrum arc still open to 
theoretical speculations, our results could be extended to support the discussion of the for- 
mation of galaxy overdense regions of a gas of massive degenerate fermions in hydrostatic 
and thermal equilibrium at finite temperatures. 
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FIG. 1: Growing mode of the density contrast 5 for DFG, radiation and CDM at super-horizon 
approximation (i. e krj 1) in dependence on the scale parameter, a. The plots correspond to 
analytical solutions in the longitudinal gauge for radiation- and matter-dominated universes. 
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FIG. 2: Analytical results for the growing mode of the DFG density contrast, 5, during the matter- 
dominated era as function of the scale parameter, a. The calculations were performed in the 
longitudinal gauge for two different scales parameterized by the wavenumber k; A: = 0.1 MPc~^ 
and A: = 1 MPc"^ 
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FIG. 3: Analytical results for the growing mode of the DFG velocity divergence, 9, during the 
matter-dominated era as function of the scale parameter, a. The calculations were performed in 
the longitudinal gauge for A; = 0.1 MPc^^ and k = 1 MPc~^. 




FIG. 4: Numerical results for the density contrast, 6, (in the synchronous gauge) for CDM, baryons, 
DFG and photons as function of the scale parameter, a. The calculations were performed in the 
longitudinal gauge for two different scales parameterized by the wavenumber k; k = 0.1 MPc^^ 
(left) and k = l MPc"^ (right). 
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FIG. 5: Comparison between the massive neutrino density contrast and the DFG density contrast 
in dependence on the scale parameter, a. It is assumed the same initial and final averaged densities 
for both fluids (the last two columns on table I). For completeness, the wavenumber here assumed 
are k = 0.1 MPc^i (left) and = 1 MPc"^ (right). 
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FIG. 6: Matter power spectrum P{k) according to the cosmic inventory described by the compo- 
nents of table I in the KCDM background scenario. 



28 



